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ABSTRACT 

In differentially rotating discs with no self-gravity, density waves cannot propagate 
around the corotation, where the wave pattern rotation speed equals the fluid rotation 
rate. Waves incident upon the corotation barrier may be super-reflected (commonly 
referred to as corotation amplifier), but the reflection can be strongly affected by 
wave absorptions at the corotation resonance/singularity. The sign of the absorption 
is related to the Rossby wave zone very near the corotation radius. We derive the 
O ' explicit expressions for the complex reflection and transmission coefficients, taking 

^ . into account wave absorption at the corotation resonance. We show that for generic 

^ ' discs, this absorption plays a much more important role than wave transmission across 

the corotation barrier. Depending on the sign of the gradient of the vortensity of the 
disc, C = K^/(2ilI]) (where VI is the rotation rate, n is the epicyclic frequency, and S 
is the surface density), the corotation resonance can either enhance or diminish the 
' super-reflectivity, and this can be understood in terms of the location of the Rossby 

wave zone relative to the corotation radius. Our results provide the explicit conditions 
CO ' (in terms of disc thickness, rotation profile and vortensity gradient) for which super- 

reflection can be achieved. Global overstable disc modes may be possible for discs with 



' super-reflection at the corotation barrier. 

, Key words: accretion, accretion discs - hydrodynamics - waves - instabilities 

o 



1 INTRODUCTION 

Differentially rotating fluid discs, ubiquitous in astrophysics, are known to exhibit rich dynamics and possible instabilities (e.g. 
Papaloizou & Lin 1995; Balbus & Hawley 1998). While local instabilities, such as Rayleigh's centrifugal instability (for discs 
with specific angular momentum decreasing outwards), gravitational instability (for self-gravitational discs with too large a 
surface density, or more precisely, Toomre Q 1), and magnetorotational instability (for discs with a sub-thermal magnetic 
field), are well understood (at least in the linear regime), global effects and instabilities are more subtle, since they involve 
couplings and feedbacks of fiuid at different locations (see Goldreich 1988 for an introduction/review). A well-known example 
is the corotation amplifier (e.g. Mark 1976; Narayan, Goldreich & Goodman 1987), which arises from the interaction across the 
corotation between waves carrying opposite signs of angular momentum. Much stronger corotation amplifications (WASER - 
wave amplification by the stimulated emission of radiation, and SWING amplifiers) can be achieved for self-gravitating discs 
(e.g., Goldreich & Lynden-Bell 1965; Julian & Toomre 1966; Lin & Lau 1975; see Shu 1992 for a review). Another well-known 
example is the Papaloizou-Pringle instability in finite accretion tori (confined between two free surfaces), in which coupling 
between waves inside the corotation with those outside, combined with refiecting inner and outer boundaries, leads to violent 
overstable modes (Papaloizou & Pringle 1984; Goldreich et al. 1986). Recent works on global disc instabilities include the 
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Rossby wave instability (for discs with a strong enough density or vortensity bump; Lovelace et al. 1999; Li et al. 2000) and 
the accretion-ejection instability (for magnetized discs; Tagger & Pellat 1999, Tagger & Varniere 2006). 

In this paper we are interested in 2D fluid discs without self-gravity and magnetic field. For disturbances of the form 
gim(/>-2t^t^ where m > and uj is the wave (angular) frequency (and thus uip — ui/m is the pattern frequency), the well-known 
WKB dispersion relation for density waves takes the form (e.g., Shu 1992) 

(lj — mO,)^ — ilP' = i<? + k^c? , (1) 

where Q. is the disc rotation frequency, uo = uj — mO. is the Doppler-shifted wave frequency, k is the radial epicyclic frequency, 
kr is the radial wavenumber, and c is the sound speed. Thus waves can propagate either inside the inner Lindblad resonance 
radius riL (defined by iD = — k) or outside the outer Lindblad resonance radius roL (defined hy Cj = k), while the region 
around the corotation radius Tc. (set hy u = 0) between riL and roL is evanescent. Since the wave inside riL has pattern speed 
Up smaller than the fluid rotation rate Q., it carries negative wave action (or angular momentum), while the wave outside roL 
carries positive wave action. As a result, a wave incident from small radii toward the corotation barrier will be super-reflected, 
(with the reflected wave having a larger amplitude than the incident wave amplitude) if it can excite a wave on the other 
side of the corotation barrier. If there exists a reflecting boundary at the inner disc radius rjn, then a global overstable mode 
partially trapped between rjn and riL will result (see, e.g. Narayan et al. 1987 for speciflc examples in the shearing sheet 
model, and Goodman & Evans, 1999 and Shu et al. 2000 for global mode analysis of singular isothermal discs). 

The simple dispersion relation Q, however, does not capture an important effect in the disc, i.e., corotation resonance 
or corotation singularity. Near corotation ^ the WKB dispersion relation for the wave is [see equation (|17|) below] 

where k^f, = m/r and S is the surface density, and the subscript "c" implies that the quantity is evaluated at r = rc. The 
quantity 

k'' _ (V X Up) ■ z 

is the vortensity of the (unperturbed) flow (where uo is the flow velocity). The dispersion relation ([2} describes Rossby waves, 
analogous to those studied in geophysics (e.g. Pedlosky 1987) 0. For k^ ^ /(? and k^ 3> fc^, we see that Rossby waves 
can propagate either outside the rotation radius Vc (when dC^/dr > 0) or inside rc (when dC^/dr < 0). In either case, we have 
fer — > 00 as r — > Tc. This inflnite wavenumber signifles wave absorption (cf. Lynden-Bell & Kalnajs 1972 in stellar dynamical 
context; Goldreich & Tremaine 1979 in the context of wave excitation in discs by a external periodic force; see also Kato 
2003, Li et. al. 2003, and Zhang & Lai 2006 for wave absorption at the corotation in 3D discs). At corotation, the wave 
pattern angular speed uj/m matches fi, and there can be efficient energy transfer between the wave and the background flow, 
analogous to Landau damping in plasma physics. Narayan et al. (1987) treated this effect as a perturbation of the shearing 
sheet model, and showed that the corotational absorption can convert neutral modes in a flnite shearing sheet into growing 
or decaying modes. Papaloizou & Pringle (1987) used a WKB method to examine the effect of wave absorption at corotation 
on the nonaxisymmetric modes in an unbound (with the outer boundary extending to infinity) cylindrical torus. 

In this paper, we derive explicit expressions for the complex refiection coefficient and transmission coefficient for waves 
incident upon the corotation barrier. We pay particular attention to the behavior of perturbations near the corotation reso- 
nance/singularity. Our general expressions include both the effects of corotation amplifier and wave absorption at corotation 
(which depends on dC^/dr). We show explicitly that depending on the sign of d(/dr, the corotation resonance/singularity can 
either enhance or diminish the super-reflectivity, and this can be understood in terms of the location of the Rossby wave zone 
relative to the corotation radius. 

Our paper is organized as follows. After presenting the general perturbation equations (section 2), we discuss the the 
wave dispersion relation and propagation diagram, and derive the local solutions for the wave equation around the Lindblad 
resonances and corotation resonance (section 3). We then construct global WKB solution for the wave equation, and derive the 
wave reflection, transmission and corotational damping coefficients in section 4. An alternative derivation of the wave damping 
coefficient is presented in section 5. Readers not interested in technical details can skip Sections 2-5 and concentrate on Section 
6, where we illustrate our results and discuss their physical interpretations. Section 6.1 contains a numerical calculation of 
the wave reflectivity across corotation and discusses the limitation of the WKB analysis. We discuss how global overstable 
modes may arise when super-reflection at the corotation is present in section 7 and conclude in section 8. 



^ A Rossby wave propagating in the Earth's atmosphere satisfies the dispersion relation oj = {2k^/k^ R){dQ.z / dO), where = fc^ + fcg, 
fis = HcosS is the projection of the rotation rate on the local surface normal vector and d is the polar angle (co-latitude). 
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2 PERTURBATION EQUATIONS 



We consider a geometrically thin gas disc and adopt cylindrical coordinate system {r,<j),z). The unperturbed disc has an 
integrated surface density E(r) and velocity uo = (0,rf2,0). The flow is assumed to be barotropic, so that the integrated 
pressure P depends only on E. Self gravity of the disc is neglected. 
The linear perturbation equations for the flow read 

d 

dl'- 
d 

dl'- 

where (5E, Su and Sh = are the (Eulerian) perturbations of surface density, velocity and enthalpy, respectively. For 

barotropic flow, 5h and 5E are related by 

2i5E 



-Su + (uo ■ V)5u + ((5u ■ V)uo = -V5h, 
-5E + V • (E Su. + uo 5E) = 0, 



(4) 
(5) 



5h ■ 



(6) 



where c is the sound speed, with — dP/dTf. 

We assume that the (p and t dependence of the perturbation are of the form 

ou, oE, oh cc e , 



(7) 



where m is a positive integer, and w is the wave (angular) frequency. We presume ix) > so that the pattern (angular) speed 
of the perturbation uip = to /m is positive (in the same direction as the flow rotation). Note that we usually assume ui is real, 
except in section 3.2 (dealing with the perturbation near corotation) where we include a small imaginary part {uj = cOr + ii^i, 
with uji > 0) to represent slowly growing disturbances. The perturbation equations ©-([S]) become 



E ^, 19,.^^ , im^^ 

—iLu — oh H — — rhrour) H Eou,* = 0, 

c' r or r 

d 

—iujSur — 2^l5uj, = — -—Sh, 
or 

—iCbSuA, + -^Sur — ~^^5h, 
212 r 

where the epicyclic frequency k is given by 



r dr 



(8) 
(9) 

(10) 

(11) 



Eliminating Sur and Su^ from equations (|8))- (|10|) . we obtain a standard second-order differential equation governing Sh 
(e.g., Goldreich & Tremaine 1979): 



d^ 
dr^ 



d 
dr 



rE/ dr 



2rnn f d , nT. 
-r- m — 
dr D 



D 



Sh = 0, 



where 



Defining 



S^D/{rE), r) = S' 



-1/2 



Sh, 



we can rewrite (|12p as a wave equation 

' d^ Dm 



2 2mQ. f d , ilT. 

m 



dr 



D 



_ gl/2 ^-1/2 

dr^ 



r? = . 



(12) 

(13) 
(14) 

(15) 



This is our basic working equation. 



3 PROPAGATION DIAGRAM AND LOCAL SOLUTIONS NEAR RESONANCES 

Consider local free wave solution of the form 

Sh oc exp 



i J kr{s)ds 



For \krr\ ^ m and away from the D = region, we find from equation ([15 



, 2 D 2mQ. f d , QT, 

kl + - + hr 1" TT 

ruj \dr U 



(16) 



(17) 
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This is the general WKB dispersion relation. Away from the region where Cj = 0, this reduces to the well-known result 
kr — —Dl<? [equation ([1])]; in the vicinity of a) = this describes local Rossby waves, with [see equations ©-((S])] 

Before studying global solutions to the wave equation (|12|) or (|15|) . it is useful to consider the special resonant locations 
in the disc. These can be recognized by investigating the singular points and turning points of the wave equation (|12p or (|15p . 
or by examining the characteristics of the dispersion relation p7|) . The special radii are 

(i) Lindblad resonances (LRs), where D = or d)^ = k^, including the outer Lindblad resonance (OLR) at iIi = k and 
the inner Lindblad resonance (ILR) eX Co = ~k. The LRs are apparent singularities of equation (|12|l or (|15p - all physical 
quantities are finite at D = 0. The LRs are turning points at which wave trains are reflected or transmitted. In the presence 
of an external force, waves are launched from LRs. 

(ii) Corotation resonance (CR), where Cj = Q. In general, the CR is a singular point of the wave equation except in the 
special case of dC^/dr = at corotation. Some physical quantities (e.g., azimuthal velocity perturbation) are divergent at 
corotation. Physically, this singularity signifies that a steady emission or absorption of wave action may occur there. 

From equation (|15|) . we define the effective potential for wave propagation by 

Veff(r) = £. + 1^ + '^(±\n—\+S^^'' — S'^'^ 
rCj \dr D ) dr^ 

= T4ff,o(r-)+AV;ff(r), (19) 

where 

, , D rn^ 2mQ. f d , \ 

AKff(0 = ^Ml^in:^ + sV2dLs-/^ (21) 
ruj dr dr'' 

Clearly, wave propagation is possible only in the region where Vcff('') < 0. Figures 1-3 depict the wave propagation diagrams 
for the cases of {d(/dr)c = 0, < and > 0, respectively. We are interested in the parameter regime c/(rS7) ^ 1 and m is of 
order unity. Note that the apparent singularity in AVcfi{r) at D — can be eliminated by analysing the wave solution around 
the LRs (see section 3.1 below). Thus we also show Vcff,o('') in Figs. 1-3. 

We now consider the behaviors of the perturbations around the LRs and CR. 



3.1 Solution Around Lindblad Resonances 



Equation (llSp has an apparent singularity at the LRs, where D ^ 0. For concreteness we will explicitly examine the outer 
Lindblad resonance (OLR); a similar solution can be found for the inner Lindblad resonance (ILR). 
In the vicinity of the OLR, equation (|15p becomes 



dr^ 



V 



dr^ k 



k^-k^l 



AmQ dk 



(22) 



rCjk dr ^ 

where k^ = ~D jf? . The last term inside (• • • ) is smaller than the second term and will be neglected. Changing the independent 
variable from r to the dimensionless integrated phase 



we have 



where 



k dr. 



k" 3 fk\ ^ 



V = Vkr 



(23) 



(24) 



(25) 



and the prime denotes differentiation with respect to z. Note that near OLR, fc^ ~ C{r — tol), with C = (— c ^dD / dr)oi^ > 



a constant, we have 



We can then express k in terms of z as 



|C1/2(^_^0^)3/2 

|Ci/2(roL-r)3/2e"3-/2 



for r > roL 
for r < roL 



k = 



-C 



1/3 



1/3 



(26) 



(27) 
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Figure 1. Wave propagation diagram in Keplerian discs: A sketcti of the effective potential 14ff,o('') (solid line) and V^fi(r) (dashed line) 
as a function of r, for the case of (dC,/dr)c = 0. Waves can propagate only in the region where Vesif) < 0, i.e., where the curves are 
below the dotted line. The three special locations arc denoted by IL (Inner Lindblad Resonance), OL (Outer Lindblad Resonance) and 
CR (Corotation Resonance). The divergence in the Vefi(r) curve around IL and OL represents an apparent singularity. 




Figure 2. Same as Fig. 1, except for the case of negative vortensity gradient, {dC,/dr)c < (or v < 0). Note the CR represents a 
singularity, and the Rossby wave zone lies inside the corotation radius. 
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Figure 3. Same as Fig. 1, except for the case of positive vortensity gradient, {d(/dr)c > (or u > 0), for which the Rossby wave zone 
lies outside the corotation radius. 



Using equation (|27p in equation (|24p we have 



V" 



36z2 



= 



Equation (|28|l has two independent solutions in terms of Bessel function (Abramowitz & Stegun 1964) 

V = V^J±2/3{z) = j^^^ [±V3Ai'{-Z) + Bi'(-Z)] , 



(28) 



(29) 



where Z = {3z/2)'^''^ , and Ai', Bi' are the derivatives of the Airy functions. Instead of using ^/zJ±2/3{z), we we can construct 
two linearly independent solutions for in a form convenient for asymptotic matching: 

for 1^1 » 1 and arg(2) = 37r/2 ^ ' 



32 



Ai'(-Z)~| 



for l^l 1 and arg(2) = 



7r\i/2 / 2 



1/6 



^) Bi'(-Z), 



sin (2 + 7r/4) for |2| 2> 1 and arg(2) = 
cxp (I2I) for \z\ ^ 1 and arg(2) = 37r/2 

where ~ indicates asymptotic expansions. This gives the connection formulae for the enthalpy perturbation at the OLE[f|: 

i ys/fc exp (- /;"^ Ifc] dr) for r < roL 



(31) 



Shi 



5h2 



/ Sjk cos ^ J^^^kdr + 7r/4^ for r ^ roL 



/57fcexp(/;°^|fc|dr) 
/S/fc sin (^/J'^^fcdr + 7r/4 



for r < roL 
for r » roL 



(32) 



(33) 



2 The usage here of 3> is used as a shorthand for the range of validity for an asymptotic expansion of a local solution. The fitting 
formulae are to be used far enough away from the resonances so that the asymptotic expansion is valid, but close enough that the local 
approximation made in (22) holds. 
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The connection formulae for ILR can be similarly derived: H 

|^/57fcexp(- r Ifcldr 



Shi 



Sh2 



/Sjk COS (/J^'^fedr + 7r/4) 



/S/fcexp (X;jfc|rfr) 
^Sjk sin {J^'^^kdr + TT/ A) 



for r >• riL 
for r <ti riL 

for r » r/L 
for r < riL 



(34) 
(35) 



3.2 Solution Around Corotation Radius 

In the vicinity of the corotation radius rc, we can rewrite (|15p dropping the m^/r^ and S terms compared to the singular 
term proportional to 1/Cj, giving 



where 



d_ j£_ 

d? "^^TiE 



dlnfl 



dlnr 



r-Rc 



Rc 



77 = 0, 



Here we have introduced a small imaginary part to the wave frequency, so that ui 
disc to a slowly increasing perturbation, we require Ui > 0. Defining 

1 



2k dr, rj - 



and recognizing that k can be treated as a constant around corotation, we have 



where e = 2krcU}i/[qL0r) and 



—lb [-- — 

dx^ V 4 X + ', 



V - — — In — - 



c d 



InC 



(36) 
(37) 

iuji. To study the response of the 

(38) 

(39) 

(40) 



^qK dr 

In equation (|39|) e > 0, consistent with the initial value problem in which the perturbation is gradually turned on starting 
from t = — oo. The parameter i/ determines the width of the Rossby wave region, Ar_R = {2c/n)\h'\. For dln^/dr ~ 1/r, we 
have |!^| ~ {H/r)c {H is the disc scale height). 

Equation (I39|l is the differential equation for the Whittaker function in complex variable z — x + ie with index 1/2 
(Abramowitz & Stegun 1964). The two linearly independent solutions convenient for the connection are 



V"- = W^,l/2(z), tp+ 



(41) 



where Tq is the Stokes multiplier defined below, and z is defined in the complex plane so that arg(z) ranges from to n. 
The particular linear combinations of Whittaker functions in (|41|) are chosen so that appropriate asymptotic expansions can 
be obtained. To obtain these asymptotic expansions and the connection formulae around the corotation, one must carefully 
consider the Stokes phenomenon, which alters the form of the asymptotic expansion of a function depending on the position 
of z in the complex plane. Since the appropriate asymptotic expansions (which connect the solution H^±„,i/2 analytically in 
different regions of the complex plane) are not readily available, we relegate the discussion of the Stokes phenomenon for the 
Whittaker function around the corotation to Appendix The resulting connection formulae are 



5h. 



5h4 



'Sjk exp k dr 

^Sjk e*'"" exp (^+ J^^kdr^ + ^S/k ■ 

'S/k exp ( + kdr) 



exp 



^" e*"" exp 



(+ /-fcdr) + (1 + ^) e-^--^ exp (- /-fcdr) 



for r 3> 
for r <^ rc ■ 

for r 2> 
for r <C f'c 



(42) 



(43) 



^ Note that for the ILR, k is real for r < riL and imaginary for r > til, while for the OLR, k is real for r > roL and imaginary for 
r < roh- 

* Note that the WKB wavenumber kr in the Rossby region ranges from oo (at r = Tc) to (at r ~ r'c it Ar^) [see equation (18)] while 

the width of the region is of order /r for most Keplerian discs. Free Rossby waves tend to be sheared away by differential rotation 
(Tagger 2001). 

^ The Stokes phenomenon is present in our Lindblad resonance analysis as well, as we see, for example, the asymptotic expansion for 
Bi' {—{3z/2)'^^^) takes the form e+^^l for a,Tg{z) = 3-k/2 but sin(2: + 7r/4) (as opposed to e'^) for arg(^) = 0. The resulting connection 
formulae are well known for the Airy functions. 
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where To and T\ are the Stokes multiphers, 



27ri 27ri e 

" V{v)T{l + uy = T{-u)V{l-uy ^^^^ 



Note that for < 1, [r(±i/)]"^ = ±!/ + 7;/^ H , where 7 = 0.5772 is the Euler constant. 



4 GLOBAL WKB SOLUTIONS AND CALCULATION OF REFLECTIVITY 

In this section we consider a wave train which approaches corotation from small radii (r ^ rm). Its propagation is impeded 
by the potential barrier between riL and roL. The incident wave is partially transmitted beyond the OLR and a reflected wave 
propagates from the ILR toward small radii. We will derive the explicit expressions for the (complex) reflection coefficient TZ 
and transmission coefficient T. 

From the dispersion relation [equation ([TJ, or equation (|f 7[) away from corotation], we find that the radial group velocity 
of the waves is 

Thus the sign of Cg/cp (where Cp = u/kr is the phase velocity) is positive for r > roL and negative for r < riL. This implies 
that in the r > roL region, the outgoing (transmitted) wave has the form exp(i kr dr) (assuming kr > Q). In the r < riL 
region, the incident wave (propagating from small radii toward corotation) has the form exp(— i kr dr), while the reflected 
wave has the form exp(i kr dr). 

A well-known property of density waves is that for r < riL the wave carries negative energy (or angular momentum), 
while for r > roL the wave carries positive energy. An incident wave exp(— i krdr), carrying energy of the amount ( — I), 
will give rise to a reflected wave TZexp{i kr dr) and a transmitted wave Texp(i kr dr). Let Dc be the energy dissipated 
at the corotation. Then energy conservation gives —I — ( — 1)|7?.|^ + |Tp + Vc, or 

i7^j^ = f + |T|^+D,. (46) 

Because of the singularity at corotation and the associated energy absorption, we first consider the simple case where the 
corotation singularity is neglected (section 4.1) before examining the general case (section 4.2). 



4.1 Neglecting Corotation Singularity: Corotation Amplifier 



Here we consider the case where the vortensity has zero slope at corotation, i.e., {d(/dr)c = 0. This would occur for disc 
models where the specific angular momentum r^Q is constant (as in the original Papaloizou & Pringle 1984 analysis), or for 
shearing sheet approximation (as in Narayan et al. 1987). In this case, there is no corotation singularity and no absorption of 
wave energy, and we always obtain super-reflection 



I7^|• 



i + iry > 1. 



(47) 



This is the essence of the corotation amplifier. 

To derive the complex TZ and T, we assume that the outgoing wave in the r > roL region (region III in Fig[T} is given 
by [see eqs. ((32)) -p3 l) ] 



Sh ■- 



where k^ = —D/(?. The connection formulae p2|l - (|33[l then give for the evanescent zone (region II in Fig. 1): 



S/k exp \ij kdr + 



Sh 



where 



2 



exp 



k\dr \ + i V S/kcxp 



\k\dr 



exp(— 611) exp 



k\dr\ + i\J S/fc exp(+Oii) exp 



\k\dr 



011 = / 1A:| dr ■ 

Using the connection formulae at ILR [eqs. (|34|) - (|35|l ]. we find that for r < riL (region I in Fig. 1) 



/■'•OL /■''OL 

:/ \k\dr = 

J riT J Tit 



\D\ 



■ dr. 



Sh: 



— ■ — e " sm 



fcdr + ) + i2\/S'/fce®" cos 



kdr + 



(48) 



(49) 



(50) 



(51) 
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Expressing this in terms of traveling waves, and defining y = JJj^ kdr — 7r/4 we fiave 



Sh~i^/S/k 
Thus the reflection coefflcient is 



7^ : 



ee„ _ ig-eii ' 



Comparing equation (|52p with equation (|48|l . we obtain the transmission coefflcient 

— i 



T : 



e+eii _ ie-eii ' 



(52) 



(53) 



(54) 



As expected, |7^p = 1 + |Tp > 1. 



4.2 Including Corotation Singularity 



As noted before, for discs with nonzero vortensity gradient (d(/dr 7^ 0), the singularity at corotation implies the absorption of 
wave energy (or angular momentum). Similar situations occur in geophysical wave systems (Dickenson 1968) and for waves in 
plasmas (Landau damping). In Appendix B, we discuss the toy problem of resonant tunneling which shares the similar energy 
absorption feature as the corotation singularity. Previous works on global modes in disc tori (e.g. Papaloizou & Pringle 1987; 
Goldreich et al. 1986; Narayan et al. 1987) suggest that the sign of dC^/dr determines whether the singularity acts to stabilize 
or destabilize a global mode. Here we derive the explicit expression for the reflectivity and the related source term T>c.- 

As in section 4.1, we assume an outgoing wave in Region III, and the connection formulae at the OLR then give for 
Region lib of Fig. [H or Fig. [3] 



5h'. 



fSjk 



exp(-enb) exp 



k\dr\ + i\/ 5//cexp(+0iib) exp 



\k\ dr ) , 



(55) 



where 



eiib= / \k\dr. (56) 
Equation (I55p is the asymptotic solution away from the in region lib. The corresponding general solution around rc is 



5h = 



ISjk 



exp(-eiib)V'+(7') + i\/ S'/fcexp(9iib)'!/'-(»-), 



(57) 



where '(/'+ and ip- are given by (I41|l . Using equations (I42p and (|43|l to match asymptotes over the corotation singularity, we 
obtain in Region Ila 



5/i 



i exp(-eiib)(l + ^roTi) + ^Ti exp(+eiib) 



iSlk e 



exp 



\k\ dr 



[exp( 



+9: 



To exp(-enb) \/ Sjk e"" exp 



\k\dr 



Using the connection formulae at the ILR, we have for Region I: 



5h 



exp(-eii) ( 1 + -ToTi ) + ^Ti exp(ei 



eiia) 



/S/k e" 



k dr 



exp(+Oii) - -To exp(eiia - 6nb] 



' b/k e cos 



kdr + 



1 + -ToTi - -Tie 



1- rp „-2e, 
-Toe 



is/sjke^- 



where y = kdr — 7r/4 and 



(58) 



(59) 



ei 



\k\dr. 



ei 



Olla + 61 



Ifcl dr. 



(60) 
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The reflection coefficient and transmission coefficient are then 



T ■ 



1 - i e-'2'^''e-2®n (1 + iToTi) - |Tie-'2'r''e-2®iia - lToe-2®iib 

l + :;(e +sin TTt/je n + — - — 

1 ^ ye -t bin niy) e ^ (r(i-i/))S 2 (r{i+i/))^ 



(61) 



(62) 



1 ^ yc ^- sm TT^^j e ^ (r(i-i/))^ 2 (r{i+i/))^ 

The dissipation term due to the corotation singularity obtained from Vc — |7?.p — 1 — |T|^ is 

ni^ cos^ TTi^ — 2€)tt ~20TTi-. I 2Tr?y ^ — 26)77^ 

-p ^ 2 (r(i+^))^g + (r(i-^))^g 

1 ^ ye -t- sm wi/j e ^ (r(i-i/))^ 2 (rci+i-))^ I 

For ^ 1 equations (|6Hl and (|62p can be simplified, and we have 



-i_ ip-©ii p+2eiib _ 1 „-2e„b 

^-^TT-^T^ + ^T -^^-^nu + 0[u% (64) 

e©n - ie-®n (ee„ _ ip-ei:^2 



peii,-e„b 

T-^^ T — 7-. 7;- 5 nv + 0[u^], (65) 

e®n-ie-©ii 2 (e©ii - ie-Qii) 

D,^^ i —4 ^ /, ,3^ Uu + 0[u% (66) 

Clearly, for ;/ = 0, equations (|6U-((65]) reduce to ((53]) and ((54)l . 



5 WAVE DAMPING AT COROTATION: ALTERNATIVE CALCULATION 

In the previous section we obtained the expression for the dissipation term Vc at the corotation resonance using the reflection 
and transmission coefficients. Here we provide a more direct derivation of T>c using the change of angular momentum flux 
across the corotation radius. 

In the absence of self-gravity, the angular momentum flux carried by the the waves in the disc is entirely due to advection. 
The time-averaged transfer rate of the z-component of angular momentum across a cylinder of radius r (in the outward 
direction) is given by (e.g. Goldreich & Tremaine 1979) 

F{r) = r^T,{r) Re[5u^(r, <;/>, t)]Re[5u^(r, <;/>, t)]. (67) 

Jo 

Using equations (f9)l- (|10|) to express Sur and 5u^ in terms of Sh, this reduces to (see Tanaka et al. 2002; Zhang & Lai 2006) 

Fir)^--^lm(sk'-^). (68) 



D \ dr ^ 

In Region III (see Figs. 1-3) the outgoing wave has the enthalpy perturbation given by (up to a proportional constant) 

5h = y/sJkT e^^ {i j ^ kdr+'^ (69) 

Calculating the angular momentum flux (setting the incoming wave flux to 1), 

i^(r > roL) =i 7^m|T|^ (70) 

we see that angular momentum is transferred outwards (positive flux) since waves in r > roL carries positive angular 
momentum. For Region I we have [up to the same proportional constant as in (|69|) ] 



5h = ^/S/k 

which gives the angular momentum fiux: 



exp ( — i / kdr + -jj + 7?.exp j kdr — ^ 



(71) 



F(r<rIL)~7^m(|7^f -1). (72) 



We see that the incident wave carries negative angular momentum outward, and the reflected wave transfers positive angular 
momentum. The net angular momentum transfer is positive (in the outward direction) for |72.| > 1. 
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Now consider the angular flux near the corotation radius, at r = (just inside corotation ) and at r = (just outside 
corotation). Integrating equation (|36p across the singularity, we flnd the discontinuity in the enthalpy perturbation derivatives: 

dSh 
dr 



dSh 


2Tvi 












r7 ' ~ 



2-nvi — 5h 
c 



(73) 



Here we have chosen to integrate from to r by going through the upper complex plane, to be consistent with the physical 
requirement of a gradually growing perturbation, turned on at t = — oo. Thus the change in the angular momentum flux 
across the corotation is 

(74) 



AFc = F(r+) - F{r, ) = \Sh'^ 



cD 

The wavefunction around rc is given by equation (|57p multiplied by T. Noting that 

1 



i>+irc) 



W.,i/2(0) 



r(i-^)' 

W_.,V2(0) + iroW.,i/2(0) = ^^^y, 



(75) 
(76) 



where we have used the convenient identity r(i^)r(l — — n/simvi^. We can evaluate the enthalpy perturbation at the 
corotation, giving 

1 _A„, COSTT!/ . A,„ 1 



Sh{rc 



'S/kT 



The change in angular momentum flux across the corotation is then 



-TrmlTI V 



n COS TTV -201 



+ le 



+ 



27r o, 



.2(r(i + /.))2 m-i^w 

With F{r < hl) = F(r-), F{r > roh) = F(r+), and thus F{r < Hl) = F{r > tol) - AF^, we find 



TV COS niy -20,- 
e 



.2(r(i + i.))2 

This expression exactly agrees with Vc given in section 4.2. 



+ 



27r 



(77) 



(78) 



(79) 



6 RESULTS AND DISCUSSION 

The key new results of this paper consist of the analytical expressions for the refiection coefficient TZ, transmission coefficient T 
and the dissipation term Dc when a wave impinges upon the corotation barrier from small radii. These expressions, (|61|l - (|66p 
and (|79|) . can be applied to discs with generic rotation and surface density profiles. 

For definiteness, here we illustrate our results using a (Newtonian) Keplerian disc model with 

J7 = Kcxr-^/', Ecxr-P, -^=(3, (80) 

ril 

where p and l3 are constants. The important parameter that determines the behavior of the corotation singularity is [see 
eq. dinj] 

Clearly the models are scale-free, and TZ, T and Vc depend only on the two parameters p and (3 (or u). Figure 4 depicts \TZ\ 
as a function of l3 for different values of p, while Figure 5 shows \TZ\, and |T| as a function of f for (3 — 0.05 and (3 = 0.1. 

A key result of paper is that wave absorption at the corotation resonance plays an important role in determining the 
reflection and transmission of waves across the corotation barrier. Without corotation resonance (as for discs with zero 
vortensity gradient, or ly — 0), super-reflection is always achieved, but \TZ\'^ — 1 ~ exp(— 20ii) (assuming Qu 3> 1, where 
©II = ©iia + 6iib = \k\ dr, with |fc| = |D|^/^/c) is rather smaU. In the presence of wave absorption at the corotation 

resonance (when 7^ 0), we find (assuming Biia ^ 1 and l^] <^ 1), 

|7^|^ - 1 ~ cxp(-2eii) + 2Tviy cxp(-2eiia), (82) 

and super-reflection can be much more prominent. From Fig. 5 we see that the transmission coefficient is generally much 
smaller than Vc. Thus the reflectivity depends mainly on the wave damping at corotation. 

Equation (|79p or (|66p clearly shows that for i/ > 0, the wave damping term due to corotation resonance is positive, Dc > 0. 
This can be understood from the fact that for i/ > the Rossby wave region lies outside r^ (see Fig. [3)l , and the dissipation 
at the corotation singularity carries away positive energy (just like the transmitted wave) so that energy conservation [see 
eq. (|46p ] requires —1 = — 17?.|^ + -I- [©d- In this case the corotation singularity enhances super-reflection, as seen in 
Figs. 4-5. On the other hand, for v < 0, the Rossby wave zone lies inside rc (Fig.[2| and the dissipation carries away negative 
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0.1 



0.05 - 



_L 



-0.05 - 



-0.1 




Figure 4. The reflection coefficient as a function of /3 = c/(rf7) for Keplerian discs with surface density proflle S ex r Note that for 
p = 1.5, wave absorption at the corotation resonance is absent. 

energy (like the reflected wave) so that —1 = — + |T|^ — \T^c\- this case the corotation singularity tends to reduce 
super-reflection, and there is a competition between the effect of the corotation amplifier [the first term in eq. (|64p ] and the 
effect of the corotational absorption [the second term]. To obtain \TZ\ > 1 we require 



u > 



1 



27J- e2©iib 



2tt 



(83) 



where the second inequality applies in the limit of Ouh ^ 1. This puts a constraint on the disc thickness and the specific 
vorticity slope (note that the sound speed c enters into both u and Quh) in order to achieve super-refiection. For example, for 
a given c, the inequality (|83|l determines the critical value of p for which \TZ\ — 1. 

Figures 4-5 also reveal an intriguing oscillatory behavior of the reflection, transmission and damping coefficients. This 
non-monotonic behavior may be qualitatively understood from the Rossby wave zone around (see Figs. 2-3). The WKB 
wavenumber kr in the near vicinity of is given by 



kj- — T" ■ 



2mQ d , , K' 
ruj ar c 



(84) 



For V > Q, the Rossby wave zone lies between Vc and Vc + Ar/j, where Atr — 2f5vrc- For quasi-normal modes to be "trapped 
in" the Rossby wave zone they must obey the the Sommerfeld "quantization" conditiorjj J^"'^^^ '^dr kr = nu ^ mv + 7r/2, 
where n = 0, 1, 2, • • • . Thus when ~ n + 1/2, the wave propagating in the Rossby zone (Fig. 3) is maximally reffected 
back to the singularity, leading to maximum negative damping and enhanced net reflection, as seen in Fig. [5] For v < Q the 
Rossby wave zone lies between — lAr^j and (Fig. 2), the wave in the Rossby zone is mostly absorbed at the corotation 
singularity. The asymmetry in the f > case and the v < Q case can also be seen in the toy problem of resonant tunneling 
(Appendix B). Note that for thin Keplerian disks, is much less than unity for reasonable density proflles (|p| ~ 1), so this 
non-monotonic behavior is of no practical interest. 



^ The phase factor 7r/2 arises from a detailed analysis of the wave behavior at r ■ 
— 37r/4 and the latter gives 7r/4. 



: rc and at r - 



Ar^: The former gives a phase of 
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Figure 5. The reflection coefficient, wave damping coefficient and transmission coefficient as a function of u for two different values of 
f} = c/(rO). Note that for u = (no corotation resonance), wave damping is zero (X'c = 0) and |7?,p — 1 is positive. 



6.1 Numerical Calculation of Reflectivity and Improved WKB Treatment 



Our analytical expressions derived in section 4-5 are based on global WKB analysis and involves several approximations. In 
particular, in our treatment of the corotation resonance (section 3.2), we assumed c <^ rQ — if this is not satisfied, some of 
the neglected terms must be included in equation (|36p and the numerical values of our solution may be modified. Thus our 
results for /3 ^ 0.1 — 0.2 should be treated with caution. 

To assess the validity of our WKB analysis, we also carry out computation of the reflectivity by numerical integration of 
equations (8)-(10). The outgoing wave boundary condition, equation (|48p . is imposed at some radius rout ^ roL, such that 

^ - (85) 



5h' 



ik + 



Sh, 



2S 2k 

where ' specifies derivatives with respect to r. At some inner radius rin <C fiL, the solution takes the form of equation (71) 
The reflection coefficient can be obtained from 



[(i-II)-H 



ik \ 5h-5h' 



m-ll)+ik] Sh-Sh' 



while the transmission is 



|T| = 



2k 



\5hl 



(86) 



(87) 



Because of the singularity at corotation, we include a small positive cui = Im(cij) for the frequency so that singularity can be 
avoided. Figure |B] depicts an example calculation for a Keplerian disc with P — 0.1 and p = 2 (so that u — 0.033). 

Figure[7]shows our numerical result compared with the calculation from the WKB analysis. There is qualitative agreement 
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Figure 6. Numerical calculation of wave reflection and transmission across corotation for a Keplerian disc with u = 0.033 and /3 = 0.1. 
The top and middle panel show the enthalpy perturbation 5h, with the solid lines depicting the real components and the dashed lines the 
imaginary components. Note that there is a discontinuity in the derivative of &h at the corotation resonance, resulting in an absorption of 
flux at the corotation. The bottom panel shows the angular momentum flux carried by the wave [see eq. (67)] . Note that F{r) is conserved 
away from corotation and the discontinuity in F at r = rc results from wave absorption. The Lindblad and corotation resonances are 
indicated by three vertical dotted lines. 



between these results. In particular, both the numerical and WKB results show that wave absorption at corotation plays the 
dominant role in determining \TZ\, and wave transmission is unimportant even for small (but nonzero) u. 

However, the WKB solution matches the numerical result closely only for \v\ <^ 1, indicating that our WKB analysis can 
be improved. For the Keplerian disk considered here, the variation in i/ was achieved by changing the background density 
index p. The WKB results shown in Figs. [5]|6] assume that the quantities Bna and Ouh are not effected by changing u. However 
p plays a non- negligible role in determining the effective wave number away from the corotation and Lindblad resonances. 
Indeed, in obtaining equation (|22p or (|36|l from p5|) . we have neglected several terms that are negligible near rc or til/ql, but 
nevertheless important away from these resonances. Noting that our connection formulae [equations (32)-(35) and (42)-(43) 
involve the asymptotic expansions of local solutions around the resonances, we can improve our WKB results by adopting the 
following ansatz: we modify the integrands in Ona and Ouh to include these dropped terms. 



e = 



,2 



D 



en ' 



4^.2 



D 



1/2 



dr 



(88) 
(89) 



Note that in the expression above, we have left out the singular term (oc ti"^) at corotation and the dominant double singular 
term (oc D~^) at the Linblad resonances since they have already been accounted for by the connection formulae in sections 
3. As shown in Fig. [SJ the improved WKB result matches the numerical solution for a much larger range of v. The increased 
values of Ona and Oub for larger \p\ act to suppress super- reflect ion, and drive reflection coefficient toward \TZ\ = 1. 
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Figure 7. The reflection coeflicient as a function of u calculated using the WKB, improved WKB, and numerical methods for a scale 
free Kcplcrian disk with sound speed c = 0.05f2r. 



7 GLOBAL OVERSTABLE MODES 

As mentioned in Section 1, global disc instabilities are often related to the super-reflection at the corotation barrier. To 
illustrate this, we consider a simple boundary condition 



5ft(nn) = 



(90) 



at the inner radius of the disc, rin <^ riL. The outgoing boundary condition at r > roL implies that the solution in the region 
r < riL is given by 

(91) 



(92) 



5h = ^/S/kexp l-i kdr + n/4 j +TZ^S/kexp I i / fcdr - 7r/4 
Letting 71 — j72.|e'^, and applying the boundary condition (|9H) yields the eigenvalue condition: 

tan (^^"^'^ kdr - n/i - <^/2^ = 
Noting that for complex eigenvalue uj = ujr + icui, the wavenumber k is also complex 

k — kr + iki — — \J (ujr + iuji)^ — ~ — \/ — + i- ^~ ^ 



where we have assumed \LOi\ <^ \iOr\. Cj,. — ui^ — mO. (< for region I) and thus the real part of the eigenvalue condition gives 

-sj{u}r — mQ.y — K.'^dr — 7r/4 — = nvr , (94) 



(95) 



where n is an integer. The imaginary part of (|92p gives (assuming \ki \ <C 



kidr ■ 



\7Z\ - 1 

|7^| + 1 
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from which we find the growth rate 



\n\ 



\Cjr\/c\/ Lu'^ — fi'-^dr 



Thus the modes are overstable {uji > 0) for \TZ\ > 1, and stable {uji < 0) for \TZ\ < 1. 



(96) 



8 CONCLUSION 

In this paper we have derived exphcit expressions for the reflection coefficient, transmission coefficient and wave absorption 
coefficient when a wave is scattered by the corotation barrier in a disc. These expressions include both the effects of corotation 
amplifier (which exists regardless of the gradient of the vortensity = K^/fiE of the background flow) and wave absorption 
at the corotation (which depends on d(/dr). They demonstrate clearly that the corotation wave absorption plays a dominant 
role in determining the reflectivity and that the sign of dC,/dr determines whether the corotation singularity enhances or 
diminishes the super-reflectivity. Our result can be understood in terms of the location of the Rossby wave zone relative to 
the corotation radius. We also carried out numerical calculations of the reflectivity. Our result provides the conditions (in 
terms of disc thickness, rotation profile and surface density profile) for which super-reflection is achieved and global overstable 
modes in discs are possible. 

In future works we will explore global oscillation modes and their stabilities in a variety of astrophysical contexts, ranging 
from accreting white dwarfs to accreting black hole systems. The possible overstabilities of these modes are directly linked to 
the effects studied in this paper and may provide explanations for some of the quasi-periodic variabilities observed in these 
systems. 



ACKNOWLEDGMENTS 

DL thanks Peter Goldreich and Ramesh Narayan for useful conversations. This work has been supported in part by NASA 
Grant NNX07AG81G, NSF grants AST 0707628, and by Chandra grant TM6-7004X (Smithsonian Astrophysical Observatory). 



APPENDIX A: STOKES PHENOMENON AND THE MATCHING CONDITION ACROSS THE 
COROTATION SINGULARITY 

As we saw in section 3.1 the perturbation equation near the corotation resonance can be solved in terms of the Whittaker 
function. However in order to provide the matching conditions we must carefully consider the effect of the Stokes Phenomenon 
on the asymptotic expansions. 

Stokes phenomenon causes the functional form of the asymptotic expansion of an entire function to be different at different 
points in the complex plane. The general asymptotic solution of the Whittaker equation for \z\ ^ 1 can be written as a linear 
combination of two functions 

P{z) = e'''^z-', Q(2) = e-'/'^z", (A-1) 

i.e., AP[z) -\- BQ{z). However, because of the Stokes phenomenon, the coefficients A and B can change when crossing the 
Stokes lines. For Whittaker functions, the Stokes lines are the positive real axis (where P is dominant and Q is sub-dominant) 
and negative real axis (where Q is dominant and P is sub-dominant). 

Consider specific solution to the Whittaker equation, with the asymptotic expansion {\z\ 3> 1) given by 

F{z) AP{z) + BQ{z), (for arg{z) = 0) (A-2) 

on the real axis. Our goal is to derive the expansion coefficients of F{z) on the negative real axis (arg(2:) = tt). To achieve 
this, we use the general results obtained by Heading (1962): 

(1) "The coefficient of the subdominant term after crossing the Stokes line = the coefficient of the subdominant term 
before crossing the Stokes line + TnX the coefficient of the dominant term on the Stokes line." 

(2) "The coefficient of the subdominant term on the Stokes line = the coefficient of the subdominant term before the 
Stokes line +^TnX the coefficient of the dominant term on the Stokes line." 

Here r„ is the Stokes multiplier for crossing the Stokes line SLig(z) — nvr in the direction of increasing arg(z) given byQ 

2j^jg-27rinsi/ 

^" = r{su)r{i + su) ^^"^^ 



Note that a typo in equation (18) of Heading (1962) has been corrected here. 
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F{z) ~ AP{z) +{B- -ToA ] Q{z), (for - tt < arg(z) < 0). 



Then in region < arg(z) < tt we have 



F{z) ~ AP{z) +[B + -ToA ) Q{z), (for < arg(z) < vr). 



Thus on the negative real axis, we obtain 



F{z 



A + l-Ti[B+ ^ToA 



P{z) +(B + -ToA ] Q{z), (for arg(2) = n). 



where s = (—1)". Note that since the pole of equation (|39p lies below the real axis, arg(2;) increases as we we move along the 
contour from z positive and real (arg(2:) — 0) to z negative and real (arg(z) = tt). 

Since equation (IA-2[) is given on one of the Stokes line, we first determine F{z) in the region — tt < arg(z) < 0: 

(A-4) 



(A-5) 
(A-6) 

(A-7) 
(A-8) 

(A-9) 
(A-10) 



The connection formulae for the two independent Whittaker functions for j2:j 3> 1 are 



Q(^) 



TiP(.) 



arg(2:) — 
arg(2:) — TT 



^4 



^W_,,i/2(2e-'") + -roW,,i/2(2). 



P(z) = e 



+ z/2-i/log{z) 



where P{z) and Q{z) have the asymptotic behavior 

e'^'exp (^J'^kdr 
for \z\ ^ 1, since in this limit j^j 3> log \z\. 



-z/2 + iylog{z) 



P{z) 

\ToQ{z) + (1 + iToTi) P{z) 



for r 3> 
for r <^ Tc 

for r Tc 
for r Tc 



a.Tg{z) = 
arg(z) = TT, 



APPENDIX B: RESONANCE TUNNELING 

To help understand the effect of wave absorption at the corotation singularity we consider the following toy problem. For a 
1-D quantum mechanical potential of the form 

V{x) = ^, (B-1) 

we have the time-independent wave equation 

^TT+^-^V = 0, (B-2) 

where we have set the Planck constant and particle mass to unity. In equation (|B-2|I . a small imaginary part ie is added to 
X, to account for the physical requirement of a growing incoming perturbation. The sign of e must be the same as the sign of 
C, so as to give us the correct physical behavior. For concreteness we will assume C > (and thus e > 0), though the same 
calculation can be performed for C < (e < 0) with similar results. 
Let y — 2ik{x + ie), where k = \PlE, we have 



, ^-- + -J^ = 0, with /. = iC7/fc,, (B-3) 

which we recognize as the Whittaker differential equation. With C > 0, y lies in the domain 7r/2 ^ arg(j/) 37r/2. Using the 
same procedure as discussed in Appendix A, we obtain the connection formulae for the general solution of equation (|B-3|I as 



^ , AP{z) + BQ{z) for arg(^) = ^/2 

^ ^ {A^TxB)P{z)^BQ{z) for arg(2) = 37r/2. ^ ' 

Consider a wave propagating from a; = — oo and impinging on the potential. The transmitted wave at a; > has the form 

1/)+ = e''=— ^{C/fc) In |2fex|+,rc/2fc a; > 1/fc. (B-5) 

The corresponding wave solution in the a; < region has the asymptotic form 

^^^^Vl-i-y-^ ^^^k:.-^(C|k)^^\1kA+Z■•^C|1k forx<-l/fc. (B-6) 
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This gives 

|7^| = 0, |T| = e""^/* (B-7) 

for a wave incident from x <Q. 

Now consider a wave incident from the a; > region towards small x. The transmitted wave is given by 

^_^^-y/2y+. ^^-ik.+i(c/m^\2k.\-3.c/2k fora;<-l/fc. (B-8) 

Connecting to the x > Q region, wc have 

^_ _^ e~^^^y'^" — T\e^^^y~" = g-»'=^+»('^/'=)ln|2fca;|-7rC/2fc _ rp^^ikx-i(C/k)\n\2kx\+irC/2k £^j. ^ ^ -^^j^ (B-9) 

This gives the reflection and transmission coefficients 

\T\ = e-^^l^, (B-11) 

for a wave incident from the positive side of the singularity [a similar result was obtained by Budden (1979) in the context of 

wave propagation in cold plasma]. 

We can define the wave absorption coefficient at the singularity as 

V = \-\Tf-\nf (B-12) 

For the forward moving incident wave (+) and the backwards moving incident wave (— ), we have 

= 1 - e-'"^/^ V- = e-2-'^/'= (i - e-2-c/'=) . (B-13) 
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